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Reaction dynamics in collisions of antiprotons on nuclei is investigated within the Lanzhou quantum molec- 
ular dynamics model. The reaction channels of elastic scattering, annihilation, charge exchange and inelastic 
collisions of antiprotons on nucleons have been included in the model. Dynamics on particle production, in 
particular pions, kaons, antikaons and hyperons, is investigated in collisions of p on '°C, Ne, *°Ca and '8!Ta 
from low to high incident momenta. It is found that the annihilations of p on nucleons are of importance to the 
dynamics of particle production in phase space. Hyperons are mainly produced via meson induced reactions on 
nucleons and strangeness exchange collisions, which lead to delayed emission in antiproton-nucleus collisions. 
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I. INTRODUCTION 


The dynamics of antiproton-nucleus collisions is a 
complex process, which is associated with the mean-field 
potentials of antinucleons and produced particles in a nu- 
clear medium and with a number of reaction channels, i.e., 
the annihilation channels, charge-exchange reaction, elastic 
and inelastic collisions. A more localized energy deposi- 
tion is able to be produced in antiproton-nucleus collisions in 
comparison to heavy-ion collisions due to the annihilations. 
Searching for cold quark-gluon plasma (QGP) with antipro- 
ton beams has been a hot topic both in experiments and in 
theory calculations over the past several decades. The large 
yields of strange particles may be produced in antiproton in- 
duced reactions, which have an advantage in comparison to 
proton-nucleus and heavy-ion collisions. The particle pro- 
duction in collisions of antiprotons on nuclei has been inves- 
tigated by using the intranuclear cascade (INC) model [1, 3? 
, 4] and the kinetic approach [5]. A number of experimental 
data was nicely explained within these approaches. A self- 
consistent description of dynamical evolutions and collisions 
of antiprotons on a nucleus with transport models is still very 
necessary, in particular the fragmentation in collisions of an- 
tiprotons on a nucleus to form hypernuclei. 

Strangeness production in antiproton induced nuclear re- 
actions has been investigated thoroughly with the Giessen 
Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model [6, 
7] and the Lanzhou quantum molecular dynamics (LQMD) 
approach [8]. The production of hypernuclei is associated 
with the reaction channels of hyperons and also hyperon- 
nucleon (HN) potential. From comparison of kinetic energy 
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or momentum spectra of hyperons to experimental data, the 
HN potential can be extracted. Also the antinucleon-nucleon 
potential is able to be constrained from particle production. 
The dynamical mechanism on strange particle production can 
be explored from the analysis of reaction channels and com- 
parison to experimental spectra. 


II. MODEL DESCRIPTION 


In the Lanzhou quantum molecular dynamics (LQMD) 
model, the dynamics of resonances (4(1232), N*(1440), 
N*(1535) etc), hyperons (A, X, 5) and mesons (x, K, n, K, p, 
w ) is described via hadron-hadron collisions, decays of reso- 
nances and mean-field potentials in a nuclear medium [9-18]. 
The evolutions of baryons (nucleons, resonances and hyper- 
ons), anti-baryons and mesons in the collisions are governed 
by Hamilton’s equations of motion which read as 


yess SS 1 
p r ap; (1) 


The Hamiltonian of baryons consists of the relativistic 
energy, the effective interaction potential and the momentum 
dependent part as follows 


Hpg = bo VP? + m; + Uint + Umom. (2) 


Here the p; and m; represent the momentum and the mass of 
the baryons. The effective interaction potential Uj; is com- 
posed of the Coulomb interaction and the local potential [13- 
18]. A Skyrme-type momentum dependent interaction has 
been used in the evaluation of the local potential energy for 
nucleons and resonances. The effect of the momentum de- 
pendence of the symmetry potential in heavy-ion collisions 
was also investigated with the isospin-dependent Boltzmann 
Uehling Uhlenbeck transport model [19-21]. 

The hyperon mean-field potential is constructed on the 
basis of the light-quark counting rule. The self-energies of 
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A and & are assumed to be two thirds of that experienced by 
nucleons. And the £ self-energy is one third of a nucleon’s. 
Thus, the in-medium dispersion relation reads 


wB(Pi Pi) = 4| (Mmg + XB)? +p? + ZY, (3) 


e.g., for hyperons as = 25N /3, ea = 257/3, 5E = XN 1/3 
and 2 7 = 2 “y /3. The antibaryon energy is computed from the 
G-parity transformation of baryon potential as 


wP p) = (mg + 58) +p? +37, (4) 


with 58 = X? and 58 = -58. 
The nuclear scalar X A and vector X X self-energies are com- 
puted from the well-known relativistic mean-field model with 


the NL3 parameter (g? „=80.82, g? „=155). The optical po- 
tential of baryon or antibaryon is derived from the in-medium 


energy as Vopt(P, p) = w(p,p) — ./p? +m. A factor £ is 
introduced in the evaluation of antinucleon optical potential 
to mimic the antiproton-nucleus scattering [22] and the real 
part of phenomenological antinucleon-nucleon optical poten- 
tial [23] as 5Y = EZ" and XN = -XN with €=0.25, which 
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leads to the optical potential V>=—164 MeV for an antinucle- 
on at the zero momentum and normal nuclear density p9=0.16 
fm~>. Fig. 1 shows a comparison of the baryon and antibary- 
on energies in a nuclear medium and the optical potentials as a 
function of baryon density. The empirical A potential extract- 
ed from hypernuclei experiments [24, 25] is well reproduced 
with the approach. The antihyperons exhibit strongly attrac- 
tive potentials in a nuclear medium. The optical potentials of 
hyperons will affect the production of hypernuclei. The val- 
ues of optical potentials at saturation density are —32 MeV, 
—16 MeV, —164 MeV, —436 MeV and —218 MeV for A(X), 
E, N, A and £, respectively. 

The evolution of mesons is also determined by the Hamil- 
tonian, which is given by 


Nm 
Hy= 2, (Vie" + o(p, p). (5) 


Here the Coulomb interaction is given by 


vee = 


rae. (6) 


j=l 


where the Ny and Nz are the total number of mesons and 
baryons, including charged resonances. The kaon and anti- 
kaon energies in the nuclear medium distinguish isospin 
effects based on the chiral Lagrangian approach in [26-30] 


1/2 
WK(P;; Pi) = [m +p? — akp? — T3¢Kp)s + (bpi + t3dxpis) | + bp; + T3dKpi3, (7) 


and 


1/2 
wgPi Pi) = [m= + p? — agp? — T3cKph + (bpi + t3dxpis) | — bxp; — T3dKpi3, (8) 


respectively. Here the bg = 3/(8f**) ~0.333 GeV fm?, the 
ax and ag, are 0.18 GeV7fm? and 0.31 GeVfm?, respective- 
ly, which results in the strength of repulsive kaon-nucleon 
(KN) potential and of attractive antikaon-nucleon KN po- 
tential with values of 27.8 MeV and —100.3 MeV at satu- 
ration baryon density for isospin symmetric matter, respec- 
tively. The t3=1 and —1 for the isospin pair KHK) and 
K°(K~), respectively. The parameters cg=0.0298 GeV?fm? 
and dx= 0.111 GeV fm? determine the isospin splitting of 
kaons in neutron-rich nuclear matter. The optical potential of 
the kaon is derived from the in-medium energy as Vop:(p, p) = 


w(p, p) — The values of mķ/mg=1.056 and 
m_/ my=0.797 at normal baryon density concluded with the 
parameters in isospin symmetric nuclear matter. The effective 
mass is used to evaluate the threshold energy for kaon and 
antikaon production, e.g., the threshold energy in the pion- 
baryon collisions V/s, = my + my. 


aP? + mz. 


Based on hadron-hadron collisions in heavy-ion reactions 
in the region of 1-2 A GeV energies [13-18], we have fur- 


ther included the annihilation channels, charge-exchange re- 
action, elastic and inelastic scattering in antinucleon-nucleon 
collisions: NN — NN, NN — BB, NN — YY and 
NN > annihilation(z,7,p,w, K,K,K*,K ,¢). Here the 
B strands for (N, A, N*) and Y(A, X, £), K(K°, K*) and K(K®, 
K`). The overline of B(Y) is its antiparticle. The cross sec- 
tions of these channels are based on the parametrization of 
experimental data [31]. The NN annihilation is described 
by a statistical model with SU(3) symmetry [32], which in- 
cludes various combinations of possible emitted mesons with 
the final state up to six particles [7]. A hard core scattering is 
assumed in two-particle collisions by the Monte Carlo proce- 
dures, in which the scattering of two particles is determined 
by a geometrical minimum distance criterion d < YoO%ror/7 
fm weighted by the Pauli blocking of the final states. Here, 
the total cross section Ctot in mb is the sum of all possible 
channels. The probability of reaching a channel in a collision 
is calculated by its contribution of the channel cross section 
to the total cross section as Penh = Och/Otor. The channel is 
randomly chosen based on the weight of the probability. 
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Fig. 1. Density dependence of the in-medium energies of nucleon, hyperon and their antiparticles in units of free mass and the optical 
potentials at the momentum of p=0 GeV/c. The empirical value of A potential extracted from hypernuclei experiments is denoted by the solid 


circle [24, 25]. 
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Fig. 2. Time evolutions of some particles produced in collisions of p on 'C at incident momentum of 1 GeV/c. 


HI. RESULT AND DISCUSSION 


The emission mechanism of particles produced in antipro- 
ton induced reactions is significant in understanding contri- 
butions of different reaction channels associated with antipro- 
tons on nucleons and secondary collisions. Fig. 2 shows the 
temporal evolutions of pions, kaons, anti-kaons and hyperons 
in the reaction p+!°C at an incident momentum of 1 GeV/c. 
It is shown that the kaons are emitted immediately after the 
annihilation of antinucleons and nucleons. The secondary 
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collisions of pions and anti-kaons on nucleons retard the e- 
missions in the reaction dynamics, i.e., rN —> A, KN > nY 
etc, which lead to the production of hyperons. At the con- 
sidered momentum below its threshold energy, e.g., the reac- 
tion NN > AA (Prhreshold= 1.439 GeV/c), hyperons are main- 
ly contributed by secondary collisions after annihilations, i.e., 
nN > KY and KN > nY. 


Phase-space distributions of particles produced in heavy- 
ion collisions were used to investigate the in-medium prop- 
erties of hadrons, in particular for strange particles [26-28]. 
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Fig. 3. The transverse momentum distributions of 2° (a), K? (b), K (c) and neutral hyperons A + X? (d) produced in p induced reactions on 


12C, Ne, *°Ca and !3'Ta at incident momentum of 4 GeV/c. 
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Fig. 4. The same as in Fig. 3, but for rapidity distributions. 


It has been shown that the mean-field potentials influence the 
spectrum structures, i.e., the kinetic energy spectra of inclu- 
sive cross sections, rapidity and transverse momentum distri- 
butions of particles etc. The phase-space structure of particle 
emission in antiproton induced reactions is also investigat- 
ed in this work. Fig. 3 shows the transverse momentum dis- 


tributions of neutral particles in collisions of antiprotons on 
different targets. The hyperon emission is coupled to the s- 
trangeness exchange reactions. The yields increase with the 
mass number of the target. Shown in Fig. 4 are the rapidity 
distributions. The target effect is pronounced in the domain 
of the antiproton-nucleon frame (y = 0.46) for 2°, K? and 
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K emissions , where a heavier target increases the produc- 
tion of pions and kaons. It leads to an opposite contribution 
for antikaon production because of the strangeness exchange 
KN —> nY. The peak of hyperon emission is close to the 
target rapidity and more sensitive to the mass . 


IV. CONCLUSION 


Dynamics on particle production in antiproton induced 
reactions, in particular pions, kaons, antikaons and hyperons, 
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